Spin projection in the shell model Monte Carlo method and the spin distribution of 

nuclear level densities 



O 
O 



> 
(N 

O 

o 
o 



o 

=5 



X 



Y. Alhassid 1 , S. Liu 1 and H. Nakada 2 
1 Center for Theoretical Physics, Sloane Physics Laboratory, 
Yale University, New Haven, Connecticut 06520, U.S.A. 
2 Department of Physics, Chiba University, Inage, Chiba 263-8522, Japan 

We introduce spin projection methods in the shell model Monte Carlo approach and apply them to 
calculate the spin distribution of level densities for iron-region nuclei using the complete (P/ + 39/2)- 
shell. We compare the calculated distributions with the spin-cutoff model and extract an energy- 
dependent moment of inertia. For even-even nuclei and at low excitation energies, we observe a 
significant suppression of the moment of inertia and odd-even staggering in the spin dependence of 
level densities. 
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The spin distribution of level densities are important 
for the calculation of statistical nuclear reaction rates 
such as those in thermal stellar reactions 0] . Knowledge 
of the spin distribution is also required for the determi- 
nation of total level densities from measured neutron or 
proton resonances @>|3j], since the latter are subjected to 
spin selection rules. 

The microscopic calculation of the spin distribution of 
level densities in the presence of correlations is a difficult 
problem. It is often assumed that the spin distribution 
follows the spin-cutoff model, obtained in the random 
coupling model of uncorrelated spins of the individual 
nucleons or excitons 0. The spin-cutoff distribution is 
determined by a single parameter, an effective moment 
of inertia. The latter is often set to its rigid-body value 
and occasionally determined empirically. 

The interacting shell model takes into account both 
shell effects and correlations and thus provides a suitable 
framework for the calculation of level densities. How- 
ever, in mid-mass and heavy nuclei, the required model 
space is many orders of magnitude larger than spaces in 
which conventional diagonalization methods can be ap- 
plied. This problem was overcome by using the shell 
model Monte Carlo (SMMC) approach 0, E| to calculate 
level densities 0, ■ SMMC level densities in the iron 
region were found to be in good agreement with experi- 
mental data without any adjustable parameters 

In the SMMC approach, thermal averages are taken 
over all possible states of a given nucleus, and thus the 
computed level densities are are those summed over all 
possible spin values. Here we introduce spin projection 
methods within the SMMC approach that enable us to 
calculate thermal observables at constant spin. We first 
discuss projection on a given spin component J z , and 
then use it to calculate spin-projected expectation values 
of scalar observables. 

We apply the method to the spin distribution of level 
densities in the iron region, and compare the results with 
the spin-cutoff model. We also extract from the spin dis- 
tributions an energy-dependent moment of inertia. Sig- 



natures of the pairing phase transition are observed in 
the energy dependence of the moment of inertia. 

We first introduce a projection on the spin component 
J z = M along a fixed z-axis. The projected partition 
function for a fixed value of M and at inverse tempera- 
ture (3 is defined by Zuifi) = T^m^~^ h where H is the 
Hamiltonian of the system and 



Tr M X = 



(aJM\X\aJM) 



(1) 



a,J>\M\ 



for an operator X. Here we assumed H to be rotation- 
ally invariant, so its eigenstates \aJM) are character- 
ized by good total angular momentum J and its mag- 
netic quantum number M with M-independent energies 
E a j. The label a distinguishes between states with the 
same spin J. The M -projected partition is then given by 



Zm{P) = J2 a j>\m\ e~P EaJ . In the following we assume 
all traces to be canonical, i.e., at fixed number Z,N of 
protons and neutrons, unless otherwise stated. 

The Monte Carlo method is based on the Hubbard- 
Stratonovich (HS) transformation e~P H = J D[a}G a U a , 
where G a is a Gaussian weight and U a is the imaginary- 
time propagator of non-interacting nucleons moving in 
auxiliary fields a. Using the HS representation, the prob- 
ability to find a state with a given spin projection M at 
temperature fi^ 1 is 
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where we have introduced the notation (X a )yy = 
J D[o\W(a)X a / J D[a]W(a), and W(a) = G a \TrU a \ is 
a positive-definite function used for the Monte Carlo sam- 
pling. $ ff = Tr[/ CT /|Tr[/ CT | in © is the Monte Carlo sign. 

In general U a is not rotationally invariant, and the 
M-projected partition TtmU^ can be calculated by J z 
projection. To this end, we use the identity 



Tr(e 



U a ) = 



E 

M=-J, 
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where J s is the maximal many-particle spin in the 
model space and tpk assumes a discrete set of val- 
ues. Using the 2J S + 1 quadrature points ifk = 
it j H fc 1 y 2 (k = —J S> ...J S ), the set of discrete func- 
tions XAf(^fe) = (2J S + i'j- 1 / 2 e l ' PkM is orthonormal, 
Y,k=-j XM(Pk)x*M>(<Pk) = hiM>- This orthogonality 



relation can be used to invert 10) 



Tr M U a 



2J S 



- jr e -^ fcM Tr 
fe=-J s 



(4) 



The trace on the r.h.s. of is a canonical trace at a 
fixed particle number A (in practice we need to project 
on both N and Z), and is calculated from the grand- 
canonical traces by a particle-number projection 



Tr 



(e^'Ua) = -^r e~ iXnA det(I + U^ fe) ) . (5) 

J N * n=l 

n,k) 



Here \n = Trn/2N S are quadrature points, and U, 
e lXn e l ' Pk ^ z XJ er is the N s x TV, matrix representing the 
many-particle propagator e lXnA e %LpkJ 'U a in the N s - 
dimensional single-particle space. In particular, e lvUa 
is a diagonal matrix with elements e ltpk7na (m a is the 
magnetic quantum number of orbital a) . The canonical 
M-projected partition Tr^ U a is calculated from Eqs. (@J 
and |J5J. 

Similarly, the canonical expectation value of an observ- 
able O at fixed M is calculated from 



(O) 
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Ttm (Qe-^) 
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where Tr m E^<t is given by (J2J , and 



Tr M (OC^) = 



2J, 



1 ^ 



e-^^Tr^e^^Uo-) . (7) 



The canonical trace on the r.h.s. of Q can be calculated 
by particle-number projection. For example, for a one- 
body operator O = X)ab( a l < ^l^) a a a f> we nn d an expres- 
sion similar to Eq. (|5J) , but each term in the sum includes 
the additional factor 



T± (aia b e i X"A e i * kJ 'U <r ) 
Tr (e lx " A e i( fi k ^-U a ) 



U, 
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where here the traces are grand-canonical. 

The spin-projected partition function at fixed to- 
tal spin J is defined by Zj{fi) = Trje- 0H = 
J2 a (aJM\e~ fJH \aJM) = J2 a e ~ PEaJ and is indepen- 
dent of M. The J-projected partition can be expressed 
as a difference of corresponding M-projected partitions 
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Eq. © holds since e~^ H is a scalar operator. Using the 
HS representation © for both M = J and M = J + 1 , 
we find 



gj(g) 
Z(p) 



Tr m = jU* 
TrU„ 



Tr» 
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where Tr mU<t are calculated as before. It is also possible 
to apply the HS transformation directly in Zj and ob- 
tain Zj(p)/Z(J3) = {{TrjU a /TiU a ) ^) w /(* a ) w . This 
relation is not equivalent to Eq. IjlOl) since U a is not rota- 
tionally invariant and TrjJ7 CT ^ Trjif = jC/ (T — TrM=j+i^cr- 
However, the calculation of Trj[/ CT requires a full spin 
projection and is considerably more time-consuming than 
the M projection required in 

To calculate the spin-projected expectation value 
(O)j = Trj {Oe-P") /Trje"^ of a scalar observable 
O (e.g., the energy), we note that Trj(Oe _/3H ) = 
Tr M=J (Oe-P H ) - Tr M=J+1 {Oe-P H ). Applying the HS 
transformation, we find 



(0)j = 



Tr M =j(OU*) 
TrU a 
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where M-projected quantities are calculated as before. 

For a good-sign interaction, U a is time-reversal invari- 
ant. Since e zipkJ * is always time-reversal invariant, so is 
e lLpkJ 'U a , and its grand-canonical trace is always posi- 
tive (since the eigenvalues of the single-particle matrix 
e'Vfc^Uff come in complex conjugate pairs). When pro- 
jected on an even number of particles, Tr (e zipkJ 'U a ) re- 
mains almost always positive. In Eq. 10} we are summing 
positive numbers with coefficients e ~ lipkM , leading to M- 
projected partition Trjvf J7 CT that can be non-positive, with 
the exception of the M = case. However, this sign prob- 
lem becomes severe only above a certain /3, and typically 
occurs at smaller values of j3 as M gets larger. We en- 
counter a similar situation for the J projection with the 
J = projection having no sign problem. In practice, the 
level density at higher spin values becomes appreciable 
only at higher excitations, and meaningful spin distribu- 
tions can be extracted except for very low excitations. 

We used the spin projection method to calculate the 
spin distribution of the partition function and level den- 
sity in the presence of correlations. In particular we cal- 
culated such spin distributions for 56 Fe (an even-even nu- 
cleus), 55 Fe (odd-even), and 60 Co (odd-odd) in the com- 
plete (pf + Og 9 / 2 )-shell, and for in the range from 
to - 2 MeV" 1 using the Hamiltonian of Ref. 7]. The 
SMMC results for Zj/Z, calculated from Eq. (JTHJ, are 
shown in Fig. For temperatures T = /3 _1 < 1.5 MeV, 
an odd-even staggering is observed in the even-even nu- 
cleus 56 Fe. In particular, Zj—q/Z is significantly en- 
hanced as T decreases. No odd-even spin staggering ef- 
fect is observed in the odd-even and odd-odd nuclei. 
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FIG. 1: SMMC spin distributions of the partition function 
ZjjZ for 55 Fe (left panels), 56 Fe (middle), and 60 Co (right). 



We calculated spin-projected thermal energies (H)m 
and (H)j as a function of (3 and used the method of 
Refs. 0, 13 to obtain the level densities pm(E x ) and 
Pj(E x ) as a function of excitation energy E x . The to- 
tal level density p(E x ) was found from (H), so we could 
determine the spin distribution pj/p at fixed values of 
the excitation energy. In Fig. |3 we show the spin dis- 
tribution of pj/p at several excitation energies for 56 Fe 
(middle panels), 55 Fe (left panels), and 60 Co (right pan- 
els). The solid squares are the SMMC results, while the 
solid lines describe fits (at fixed E x ) to the spin-cutoff 
model 



j(j+i> 



(12) 



with an energy-dependent spin-cutoff parameter a as the 
only fit parameter. The spin-projected density pj(E x ) 
in l|12f> is normalized such that ^Zj(2J + ^)pj{E x ) « 
p(E x ). Equation l|12l) follows from the random coupling 
model, in which the distribution of the total spin vector 
J is Gaussian 4]. At intermediate and high excitation 
energies the spin-cutoff model seems to work well for all 
three nuclei. However, for the even-even nucleus 56 Fe, 
we observe an odd-even (spin) staggering effect below 
E x ~ 8 MeV that cannot be explained by the spin-cutoff 
model. 

The energy-dependent spin-cutoff parameter a 2 (E x ), 
obtained by fitting pj/p to Eq. 1)12(1 . is shown (solid 
squares) versus E x in the top panels of Fig. for 55 Fe, 
56 Fe, and 60 Co. The quantity a 2 can also be obtained 
from fits to pm/p [in the spin-cutoff model pu/p = 
(27TCT 2 ) _1 / 2 e _M l 2a ] but the results are similar. Despite 
the deviation from JEJl at E x < 8 MeV in 56 Fe, the fitted 
a 2 (E x ) represents well the average behavior of p.j/p. 

There is not much data available regarding the spin- 



FIG. 2: Spin distribution of level densities, pj/p, at constant 
excitation energy E x for 55 Fe, 56 Fe and 60 Co. The SMMC re- 
sults (solid squares) are compared with the spin-cutoff model 
Pi with a 2 fitted to the SMMC results (solid lines) , and with 
a calculated from Eq. 11311 using the rigid-body moment of 
inertia (dashed lines). 



cutoff parameter. A few available experimental data 
points are shown for 55 Fe The SMMC calculations 
are in general agreement with the experimental data. 
a 2 is related to an effective moment of inertia I through 



IT 



(13) 



where T is the temperature. Using Eq. (|13fl we can 
convert the SMMC values of a 2 to an energy-dependent 
moment of inertia I(E X ). The results are shown in the 
middle panels of Fig. [21 (solid squares). For comparison 
we also show the rigid-body value I/h 2 = 0.0137A 5 / 3 
MeV -1 (dashed lines), and (for 56 Fe only) half the rigid 
body value (long dashed lines) of the moment of inertia. 
In all three nuclei, I(E X ) is a monotonically increasing 
function of E x and is close to the rigid-body value at in- 
termediate and high excitations. However, for energies 
below ~ 8 — 10 MeV we observe a suppression that is 
particularly strong for the even-even nucleus 56 Fe. 

The energy dependence of the moment of inertia ex- 
tracted from the spin distributions originates in pair- 
ing correlations. To demonstrate that we calculated the 
average number of J — nucleon pairs (At A), where 

A 1 = E ama > (-l) ja - m «alrnA-rn a - The SMMC re- 
suits for proton (p-p), neutron (n-n) and proton- neutron 
(p-n) pairs are shown versus E x in the bottom panels 
of Fig. |3 The rapid decrease of the number of p-p and 
n-n pairs for 56 Fe is strongly correlated with the rapid 
increase observed of the moment of inertia. The correla- 
tion between / and the number of pairs suggests that 
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FIG. 3: Shown from top to bottom are the spin-cutoff param- 
eter a 2 (extracted from the spin distributions), the moment 
of inertia I in il l Ml and the number of J = pairs (A^A) 
for 55 Fe (left panels), 56 Fe (middle), and 60 Co (right). The 
SMMC results for a 2 and I are denoted by solid squares. The 
open circles with horizontal errors in the a 2 panel of 55 Fe are 
experimental data 9] . The dashed lines correspond to a rigid- 
body moment of inertia while the long dashed lines (for 56 Fe 
only) correspond to half the rigid-body moment of inertia. 

at low excitation energy the nucleons behave as con- 
densed BCS pairs, leading to moment of inertia values 
that are significantly smaller than the rigid-body value. 
Thus the spin distributions provide a thermal signature 
of the pairing phase transition. Thermal signatures of 
pairing correlations were previously observed in the heat 
capacity [lfj, 111) . The strong odd-even effect seen in the 
moment of inertia can be explained by a simple pairing 
model plus a number parity projection 

The strong suppression of the moment of inertia at 
low excitations for even-even nuclei has a signature in 
the J = level density. In Fig. 21 we show the total and 
J = level densities for 56 Fe. The SMMC total density 
(solid squares) is in good agreement with the backshiftcd 
Bethe formula (BBF) (dotted-dashed line). The SMMC 
J = level density is shown by the solid circles. The solid 
line describes the J = level density obtained from l|12|) 
(and the BBF for the total p) using the energy-dependent 
moment of inertia in Fig. 01 T ne dashed line is also from 
Eq. i|12[) but with a rigid-body moment of inertia. The 
rigid-body curve agrees with the SMMC results at high 
excitations, but shows deviations for energies below ~ 
8 MeV. The enhancement of the J = level density 
for E x < 8 MeV reflects the decrease of the moment of 
inertia due to pairing correlations. 

In conclusion, we have introduced spin projection 
methods in the shell model Monte Carlo approach and 



used them to calculate the spin distributions of level 
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FIG. 4: Total and J = level densities for 56 Fe. The total 
SMMC level density (solid squares) is well described by the 
BBF level density (dotted-dashed line). The J = SMMC 
level density (solid circles) is compared to the J = level 
density inferred from (1 1 2t with the fitted moment of iner- 
tia (solid line) and the rigid-body moment of inertia (dashed 
line). The suppression of the moment of inertia at low exci- 
tations enhances the J — level density E x < 8 MeV. 

densities. The energy-dependent moment of inertia ex- 
tracted from these distributions displays an odd-even ef- 
fect that is a signature of the pairing correlations. 
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